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Abstract. - We propose a theoretical description of the sin {2ip) current-phase relation in 
SFS junctions at the O-tt cross-over obtained in recent experiments by Sellier et al. [1] where 
it was suggested that a strong decoherence in the magnetic alloy can explain the magnitude of 
the residual supercurrent at the O-vr cross-over. To describe the interplay between decoherence 
and elastic scattering in the ferromagnet we use an analogy with crossed Andreev reflection in 
the presence of disorder. The supercurrent as a function of the length R of the ferromagnet 
decays exponentially over a length ^, larger than the elastic scattering length Id in the absence 
p^ , of decoherence, and smaller than the coherence length l,p in the absence of elastic scattering 

Q_) ' on impurities. The best fit leads to ^ ~ ^j^ ' ' /3, where (,j^ is exchange length of the 

diffusive system without decoherence (also equal to ^ in the absence of decoherence). The fit 
of experiments works well for the amplitude of both the sin ip and sin {2ip) harmonics. 
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\^^ , Introduction. - When Cooper pairs from a superconductor (S) penetrate a ferromagnet 

f^ ' (F) at a SF interface the spin-up electron decreases its spin energy because of Zeeman split- 

■^ i ting whereas the spin-down electron increases its spin energy. As a consequence the spin-up 

electron accelerates and the spin-down electron decelerates so that Cooper pairs acquire a 

finite center of mass momentum Afc. This induces spatial oscillations of the pair amplitude 

in the ferromagnet. If the length R of the ferromagnet is well chosen the Josephson relation 

of a SFS junction can be inverted, giving rise to a 7r-junction [2-5]. The determination of 

'"^ ' the current-phase relation in SNS (where N is a normal metal) and SFS junctions is a long 

^ I standing problem (see the recent review by Golubov et al. [6]). The oscillations of the transi- 

O ■ tion temperature of FS superlattices [7] as a function of the thickness of the ferromagnet are 

another consequence of the 7r-coupling. 

SFS TT-junction have been probed recently in experiments in which the ferromagnet is 

k>( [ a magnetic alloy with a sufficiently small exchange field [8-12] so that the period of the 

5^ I spatial oscillations of the pair amplitude is large enough. In the following we discuss a recent 

5r ! experiment [1] in which half-integer Shapiro steps were observed in a Nb/CuNi/Nb junction 

at the O-vr cross-over, indicating a sin (2ip) current-phase relation. For highly transparent 
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interfaces and without decoherence the current-phase relation is related to the derivative of 
the energy of the Andreev bound states with respect to the superconducting phase difference 
between the two superconductors [3-6]. The exchange field in the ferromagnet generates a 
Zeeman splitting of the spectrum of Andreev bound states, which was the explanation of the 
sm{2(p) harmonics at the O-tt cross-over proposed by Sellier et al. [I]. To explain the tiny 
magnitude of the residual supercurrent at the O-tt cross-over, Sellier et al. [1] suggested the 
existence of a strong decoherence in the magnetic alloy. However within this assumption the 
Andreev bound states are broadened so that the supercurrent cannot be anymore evaluated 
as the derivative of the bound state energies. This calls for a specific modeling of n junction 
involving decoherence in the ferromagnet. 

The origin of decoherence in a magnetic alloy is not well characterized experimentally, 
because it is not possible to carry out weak localization and universal conductance fluctuations 
as a function of an applied magnetic field. Quantum coherence in a ferromagnet was however 
studied in a recent experiment using time-dependent universal conductance fluctuations [13]. 
One source of decoherence in a ferromagnet is spin waves but other effects can play a role 
such as spatial heterogeneities of the exchange field or domain wall motion [14, 15]. 

We discuss here the effect of decoherence in the magnetic alloy, motivated by the exper- 
imental observations of a very small residual supercurrent at the O-tt cross-over and more 
generally to a huge reduction of the supercurrent compared to a model not involving decoher- 
ence [1, 12]. We show that decoherence in the ferromagnet implies the existence of a length ^, 
intermediate between (i) the elastic mean free path l^i due to elastic scattering on impurities 
in the absence of decoherence, and (ii) the coherence length l^p due to decoherence in the 
magnetic alloy in the absence of elastic scattering on impurities. The sin ip harmonics due 
to the coherent transfer of a charge-(2e) is proportional to exp (— i?/^), and the sin {2(p) har- 
monics due to the coherence transfer of a charge-(4e) is proportional to exp (— 2_R/^), where 
R is the length of the ferromagnet. The length scale ^ can be smaller than R even though Icp 
is larger than R. Therefore we base our description on the first two terms of an expansion 
in the coherent transfer of multiples of a charge-(2e). Similarly to Ref. [16] the dressing by 
multiple local Andreev reflections is described non perturbatively so that our description is 
valid for highly transparent interfaces. We also treat rigorously the geometrical effects related 
to propagation parallel to the interfaces. 

The models. - We consider a model in which a 3D ferromagnet is connected to two 3D 
superconductors (see Fig. n}(a)). The superconducting electrodes are described by the BCS 
Hamiltonian 

^BCS = X! ~^ (4,<tC/3,^ + cl;„Ca,a) + A ^ (4,T"=a,i + ^ajC^J j , (1) 

{a.j3),a a 

where a and /3 correspond to neighboring sites on a cubic lattice. The ferromagnetic electrode 
is described by the Stoner model 

(q,/3),(T a 

where ft-ox is the exchange field. We add to a term describing elastic scattering by impu- 
rities. The ferromagnet is connected to the superconducting layers by the Hamiltonian 
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X! ^*a(6) [^l.a,a{bf°''''^^ + (^tM,S^aM,a{b)) ' C^) 
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Fig. 1 - Schematic representation of the 3D SFS junction (a); the diagram contributing to the simp 
harmonics (b); the three diagrams contributing to the sin2i/9 harmonics (c), (d), (e). For clarity, 
renormalization by local process is not shown on these diagrams. The sites at the interfaces are 
labeled by a^ and a^ (interface a-a) and h, Pi (interface 6-/3). 



where the label a runs over all sites at the interface. 

To describe propagation in the ferromagnet we use the advanced Green's function of a 
bulk ballistic 3D ferromagnet 



1,1,A 
9a.b = 
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where l^p is the phase coherence length. A similar expression is obtained for g^'j^' . We will 
discuss how to incorporate disorder in the ferromagnetic electrode. The Green's functions of 
a superconductor can be found in the literature [17]. 

The Nambu representation of the tunnel amplitudes is given by 



f — f* — f 









-iip/4 



(5) 



and similar equations are obtained for tb^p and i^gj,. 
Supercurrent. - The supercurrent is given by 



Is^^ 



J °^ Tr {a' [L.a {Ci 



G„ 



k^{Gia-G^,a)]}dio, (6) 



where the trace corresponds to a summation over the "11" and "22" matrix elements in the 
Nambu representation, a summation over the two spin orientations (formally equivalent to 
summing over the "33" and "44" components in the spin (g) Nambu representation [18]), and 
a summation over the channel labels. Eq. Q can be derived from Keldysh formalism [19,20] 
by noting that in equilibrium the Keldysh Green's function takes the simple form G^'" = 



nF{uj) 



G^ ~G^ 



equation G — g 



. The fully dressed Green's functions Gij are determined through the Dyson 
g (g) S (g) G, where E is the self energy corresponding to the bonds i^ q, and 
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tb^p and denotes a summation over the sites a^, a^, bi, Pi (see Fig. ^ and a convolution 
over the time arguments that becomes a simple product once the Dyson equation is Fourier 
transformed to the energy variable lu. The supercurrent takes the simpler form 



^-4 






(7) 



where [] is an anticommutator and a^ one of the Pauli matrices. 



Perturbative expansion of the supercurrent. - Now we describe a perturbative expansion 
in which we include the coherent transfer of a charge-(2e) and charge-(4e) object while keeping 
a non perturbative description of local processes. The Green's fmiction Ga,a can be expanded 
according to Ga,a = J2n G'o"o, where Gi^i describes the sin (mp) harmonics due to the coherent 
transfer of a charge-(2ne) object. We obtain Ga,a = Ka^a9a.a, 



G(i) 
G(2) 



Ka,aXaMKb,bXb,a 



9b,a 



with Xn 



Ka.aXa,bKb,bXb,aKa^aXa.bKb,b 
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9bM + 



(8) 
(9) 



ya,a'^a,a9a,a^a,a-] ^a.a 



i~x„ 



Kafi, Kb^a- The channels labels are kept implicit. 



The simp harmonics of the supercurrent is given by Ig sin 99, with 



'^b,a^a^a9o,^a 

, and with similar expressions for Kb^b, 
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Tr 



ga.,bVb,bgb,aWa,a 



Tr 



ga.bVb,b9b,aW'^,<\ } , (10) 



where the sum over a and 6 is a sum over all channels at the two interfaces. The vertices V ^ 
W and W^ containing information about the dressing by local processes are given by 



a. a — ^a,a9oi,ct^Oi,a^a,a9a,,a,^a,ct [9a, a-) ^ \ ^OL^a^a^a- 



(11) 

(12) 
(13) 



The same perturbative expansion can be applied to the amplitude Ig of the sin (2i^) harmon- 
ics. The three diagrams on Fig. QKc), (d), (e) give rise to 24 terms that can all be evaluated 
explicitely. One of these terms is 



+ OC 



ai,&i,a2,f)2 



1,1 T/1,2 
9ai,bi ^fci,6i 



2,2 



ta2,1 1,1 T^l,2 2,2 T/r/2,1 

''^01,02 "02,62 ft2, 62^62, 02 12,01 ■ 



(14) 



The quantities Vj, ' ^ and / V^'^ {R)W^'^ {R)dR (with R the vector between oi and 02) are eval- 
uated by a Fourier transform. Vf,'^ is proportional to / y^^^(fc||)(ik|| a.nd J V^'^{R)W^'^{R)d'R 
is proportional to t/^'^(fc|| — 0)W'^''^{k\^ = 0). We now evaluate the average over disorder of 

9a b 9b a ^^'^ show that, like for crossed Andreev reflection [21,22], disorder changes the 
value of the exponent in the geometrical prefactor compared to the ballistic case and therefore 
enhances the supercurrent through the ferromagnet compared to the ballistic case. We also 
calculate the coherence length ^ and the effective Fermi wave-vector mismatch AK in the 
presence of both disorder and decoherence. 
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Effect of disorder in the ferromagnet. - The average Green's function as a function 

of separation d decays like g^^^' (d) — g^'i^' (d) exp {—d/ld) [17] where g^'^' (d) is given by 
Eq. Q . The mean free path Id due to elastic scattering on impurities is supposed to be much 
smaller than the ballistic coherence length l^p . The calculation of the diffusion probability 



tp- 



V{d) = J:^gl^{d)glM, 



(15) 



with ao the lattice spacing, is formally analogous to the calculation of the subgap conductance 

of a disordered superconductor [21,22], except for different phase factors. Following Ref. [21] 

we obtain 

1 / d' 



V{d) ^ - 



47r?ii?d ^""P I ^ 



with D ~ vpld/i the diffusion constant, and with 




(16) 



(17) 



We deduce from Eq. H17|l the following expression of the effective wave-vector mismatch AiC 
and coherence length ^: 



(18) 



2 
Integrating over all channels at the two interfaces we obtain 



^^^9]^M^{d)gl'l''^{d) 
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= -N, 



ao 



ch 
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tp^hD 



1/^-iAK 



(21) 



with d = y/W 



Fit of the experiments. - Now we use the model discussed above to fit quantitatively the 
recent experiments by Sellier et al. [1] . The value of the bulk hopping amplitudes tp and ts 
are chosen equal to avoid multiplying the number of parameters, and such that the diffusion 
constant of the ferromagnet is close to -D = 4 cm^s""'^ [1]. We use ts = tp = 'b x 10^ K and 
the elastic mean free path is Z^ = 10 A [12]. The Fermi wave- vector kp in the superconductor 
is chosen equal to kp — 1 A~^, and the lattice parameter is chosen equal to oq = 1 A. Since 
the interfaces are highly transparent we choose ta — ti, = ts = tp. We fix the ratio between 



Im and ^ 



(ball) 



hvp/hcx in the ballistic system to be smaller than unity, and determine the 



exchange field in such a way that the sin ip harmonics vanishes for the same value of R as in the 



experiment. We first tried a fit with L 



^(ball) 



7150 A, ft-ex = 70 K, and with ^ == 100 A. 



This fit is not compatible with the three experimental points with the largest values of R 
and the residual value of the supercurrent at the O-tt cross-over is far too large (curve (a) 
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Fig. 2 - Variation of the amplitudes 7g and Ig of the sin Lp and sin (2(p) harmonics as a function of 
the thickness R of the ferromagnet. The experimental points [1, 12] are indicated on the left panel. 
We use D = 3.5 cm^s"\ Id = 10 A. (a) corresponds to /icx = 70 K, l^ = 7150 k, ^ ^ 100 A. (b) 
corresponds to /lox = 190 K, l^ = 1300 A, .^ = 45 A. (c) corresponds to /lox = 460 K, l^ = 270 A, 
5 = 20 A. 



on Fig. [21). To obtain a better agreement with experiments we increase /lex and decrease l^ 
The fit (b) corresponds to he 



190 K, l^ = 1300 A, C 

270 A, ^ 



45 A. The fits (c) on Fig. |21 is 
even better, with /lex = 460 K, l^ = 270 A, ^ = 20 A. The values of ^ can be compared to 
the values of Q ' — y/hDjh^ for the diffusive system in the absence of decoherence. The 
values of ^ ' /^ are approximately equal to 1.5 (fit (a)), 2 (fit (b)) and 3 (fit (c)), whereas 
Q ' Y^ = 1 in the absence of decoherence. 

The amplitude of the sin {2ip) harmonics at the O-tt cross-over is 4 /xA in experiments. 
The fits (a) and (b) lead to much larger values whereas the fit (c) has the correct order of 
magnitude (see Fig. |21l. The agreement of the fits of the amplitudes of both the sin ip and 
sm(2ip) harmonics suggest the validity of the fit (c). The renormalization of the coherence 
length ^ compared to its value f j^ ' Mn the absence of decoherence may depend on the quality 
of the magnetic layer since it is expected qualitatively that decoherence is reduced with less 
inhomogcneities in the exchange field. 

Conclusion. - To conclude we have provided a modeling of the sin (2(^) current-phase 
relation in SFS junctions at the O-tt cross-over. This model is motivated by the fact that the 
residual supercurrent at the O-tt cross-over is very small, therefore suggesting that decoherence 
in the ferromagnet plays a relevant role, as suggested by Sellier et al [1]. Decoherence was 
introduced through a phenomenological coherence length l^p. Elucidating precisely its micro- 
scopic origin is difficult experimentally but several factors (spin waves, inhomogeneities in the 
exchange field, motion of domain walls) may play a role. The supercurrent was calculated 
through an expansion in the number of non local processes connecting the two interfaces. Low- 
est order processes correspond to a coherent transfer of a charge- (2e) contributing to the simp 
harmonics and the next order corresponds to a coherent transfer of a charge- (4e) contributing 
to the sin {2(f) harmonics. Like for crossed Andreev reflection [16] these non local processes 
are dressed by multiple local Andreev reflections that were included in a non perturbative 
fashion. We also included the geometrical effect of propagation parallel to the interfaces. 

We found that in the ferromagnet disorder effects (characterized by the elastic scattering 
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length Id in the absence of decoherence) and decoherence (characterized by the coherence 
length l^ in the ballistic limit) generate a new coherence length ^, intermediate between l^ and 

lip , and smaller than the exchange length ^j^ ' evaluated in the diffusive limit in the absence 
of decoherence. The amplitude of the sirup harmonics is proportional to exp (—R/^) and the 
amplitude of the sm{2ip) harmonics is proportional to exp(— 2i?/^) so that exp(— i?/^) is a 
small parameter in the expansion in non local Andreev processes. The approach is consistent 
in the sense that the fit to experiments shows that exp (— i?/^) is a very small parameter [we 
have R/^ = 8.75 for the fit (c) on Fig. [2] so that the sin (3(/?) and higher harmonics can indeed 
be neglected. 
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